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Atomic and molecular gas in the starburst galaxy NGC4945 
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Abstract. Spatial and kinematical correlations between the Hi and 12 CO(2— 1) emission of the southern spiral 
galaxy NGC 4945 are studied with a common angular resolution of ~ 23" (corresponding to 750 pc at D — 6.7 Mpc) 
and a velocity resolution of ~7kms _1 . The 1.4 GHz continuum emission is also observed. The Hi kinematics yield 
a galaxy mass of ~ 1.4 x 10 M© within radius _R=380", with molecular and neutral atomic gas each contributing 
~2%, respectively. A central continuum source of size 7 '.'6 x 3'.' 4 (250 x 110 pc) is enveloped by a molecular cloud 
of mass 1.5 x 10 9 Mq for R < 7'.'5, and is rapidly rotating with V la t ~ 160 kms -1 . Hi emission from the central 
region at velocities \V— V syB \ > 200 kms -1 may be related to optically detected gas that is believed to trace an 
outflow directed towards the halo. Nuclear Hi absorption at V - V Bys ~ +80 kms -1 suggests inflow towards the 
centre, that was so far only seen in molecular lines. Hi features at each end of the major axis (\R\ ~ 600") are 
interpreted as spiral arms that are viewed tangentially and that also cause prominent emission features in the 
radio continuum, Hi, and CO further inside the galaxy. A central elongated region showing non-circular motions 
is interpreted as a bar which fuels the nuclear starburst. The Hi and CO position- velocity data have been analysed 
using linear resonance theory, and possible locations of resonances are identified. 
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1. Introduction 

Studies of the formation and evolution of spiral galaxies 
require knowledge of their morphological and kinematical 
properties. Atomic and molecular hydrogen (Hi and H2) 
are the main interstellar gas components of spiral galax- 
ies and provide excellent tracers to elucidate spiral struc- 
ture and rotation. Since it is difficult to directly observe 
H2, carbon-monoxide (CO) is commonly used to probe H2 
column densities and molecular masses. Whereas Hi is a 
good tracer for the outer regions of spiral galaxies, the 
inner regions are often better studied in CO. 

The starburst galaxy NGC 4945 is particularly suited 
for high resolution and high sensitivity studies. Being 
a member of the Centaurus group, the nearby edge-on 
galaxy contains a Seyfert 2 nucleus and is classified as 
SB(s)cd or SAB(s)cd (de Vaucouleurs 1964; Braatz et 
al. 1997). Distance estimates range from 3.8 Mpc (de 
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Vaucouleurs 1964; Bergman et al. 1992) to 8.1 Mpc (Baan 
1985). In order to be consistent with other relevant stud- 
ies, throughout the paper 6.7 Mpc is used. 

Although at optical wavelengths the active nucleus 
is obscured (JVh ~ 10 24 cm~ 2 ; Guainazzi et al. 2000; 
Madejski et al. 2000), there is evidence for a nuclear su- 
perwind (Chen & Huang 1997; Lipari et al. 1997) and, at 
A = 100 fixa, the central region is one of the three brightest 
IRAS point sources beyond the Magellanic Clouds (IRAS 
1989). The nucleus is well-defined at radio frequencies 
(Si.4GHz = 4.9 Jy; Elmouttie et al. 1997) and shows lit- 
tle ultra compact structure on a milliarcsec scale (Preston 
et al. 1985; Sadler et al. 1995). Microwave transitions of 
molecules have been detected in dense clouds that enve- 
lope the nucleus; this includes the first discovered 'mega- 
maser' (Dos Santos & Lepine 1979; Batchelor et al. 1982; 
Whiteoak & Gardner 1986), a possible circumnuclear disk 
of radius 0.3 pc with a binding mass of ^10 6 M© (Greenhill 
et al. 1997), and numerous transitions at cm and mm 
wavelengths (see e.g. Whiteoak 1986; Henkel et al. 1990, 
1994; Curran et al. 2001). CO studies suggest the pres- 
ence of a nuclear gas ring (Whiteoak et al. 1990; Bergman 
et al. 1992; Dahlem et al. 1993). Properties of NGC 4945 
relevant to this paper are summarized in Table |l|. 
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Table 1. Properties of NGC4945 (9 mag, SB(s)cd) 



Position 


a(2000) 


13 h 05 m 27f4 




of nucleus 


(5(2000) 


-49° 28' 05" 


1, 2 


Distance 3 ^ 


D 


6.7 Mpc 


3 


Diameter 


D25m 


17' 


4 


Systemic 


vw 


563 and 561±4kms _1 


5, 7 


velocity 1 ^ 


Vlsr 


555 kms -1 


6 


Inclination 


i 


78 ±3° 


6, 7 


Position angle 


PA 


43 ± 1° ; 45 ± 2° 


6, 7 


Turn-over radius 


1 W max 


5.' 45 ; 6.'6 


6, 7 


Turn-over- velocity 


Vmax 


165 and 183 ± 4km s" 1 


6, 7 


Total mass 


Mt 


8.8 and 15 x 1O 1O M 


6, 7 



a) The value of 6.7 Mpc was chosen for consistency with 
previous Hi and CO studies; at this distance 23" corresponds 
to 750 pc. Evidence for D = 4 Mpc has been summarized by 
Bergman et al. (1992). If this value were used, gas masses and 
luminosities have to be multiplied by a factor of 0.36 except 
for dynamical masses, for which a factor of 0.6 is appropriate. 

b) Vhei = Vim + 4.7 km s^ 1 
References: 

1) Whiteoak & Bunton (1985) 5) Whiteoak et al. (1990) 

2) Elmouttie et al. (1997) 6) Abies et al. (1987) 

3) Whiteoak & Gardner (1977) 7) Dahlem et al. (1993) 

4) de Vaucouleurs et al. (1991) 

Early optical spectroscopic observations suggest that 
the rotational velocity increases linearly out to a radius of 
r = 360" but also show some evidence of non-circular mo- 
tions (Peterson 1980; Carranza & Agiiero 1983). Studies 
of the disk radio continuum, Hi, and CO emission have 
been made by Elmouttie et al. (1997), Abies et al. (1987), 
and Dahlem et al. (1993) respectively. The large-scale ra- 
dio emission has a steep (frequency) spectral index (cv ~ 
-l.lztO.l); no emission was detected from the halo. Hi was 
found to extend out to r = 7' with an additional feature at 
10' in the south-west. CO(l— 0) emission is concentrated 
towards the nucleus, but was also detected out to 6' along 
the major axis. 

We have used the Compact Array of the Australia 
Telescope National Facility (ATCA, see Frater & Brooks 
1992) to study the 1.4 GHz continuum and Hi emission 
of NGC 4945 in more detail and with higher sensitivity 
than before. We have also observed CO(2— 1) emission us- 
ing the Swedish-ESO Submillimetre Telescope (SEST; see 
Booth et al. 1989). In addition to the results presented by 
Dahlem et al. (1993) our CO(2— 1) data are not confined 
to the nuclear region but also include significant parts of 
the disk. H 1 and CO have been compared at a common 
angular resolution of ~23" and channel spacings of 6.6 and 
6.3 km s -1 , corresponding to velocity resolutions of ~7.25 
and 6.5 km s" 1 , respectively; 23" corresponds to 750 pc at 
D = 6.7Mpc (see also Table @, footnote 'a'). 

2. Observations and data reduction 



2.1. 1.4 GHz continuum emission and H 1 synthesis 
data 

The ATCA was equipped with 20 cm FET receivers with 
system temperatures between 33 and 55 K. Four observ- 
ing runs provided 53 independent baselines, and were 
sufficient to achieve a satisfactory sampling of the u- 
v plane. Table || displays the configurations, observ- 
ing dates and baseline ranges. A standard correlator 
setup of 512 channels was used with a total bandwidth 
of 8 MHz centred on 1418 MHz. The channel separa- 
tion was 15.6 kHz (3.3 km s -1 ). An observing cycle was 
adopted in which a 40-minute tracking of NGC 4945 was 
bracketed by alternate 5-min observations of the phase- 
calibration sources PKS 1215-457 and PKS 1320-446. The 
radio galaxy PKS 1934-638 was observed as a flux-density 
calibrator (assumed to have a flux density of 16.4 Jy at 
1.4 GHz). Observations of the quasar PKS 0407-658 pro- 
vided antenna spectral bandpasses. 

Table 2. ATCA configurations for the Hi observations. 



Config. 


Date 


Baselines 


[m] 


1.5B 


Oct. 18, 1991 


31 


to 


4301 


1.5C 


Dec. 3, 1991 


77 


to 


4500 


0.75 C 


Aug. 22, 1992 


46 


to 


5020 


6.0C 


Feb. 12, 1993 


153 


to 


6000 



The data were reduced with an ATNF modification 
(Killcen 1992) of the Astronomical Image Processing 
System (AIPS) of the US National Radio Astronomy 
Observatory. A single data file of correlated u-v data was 
produced for each observing run. After editing of the 
spectra obtained with individual baselines for phase and 
amplitude errors, the central 400 spectral channels were 
selected. Those channels free of Hi contributions were 
averaged in the u-v domain for each of the four data 
sets and merged into a continuum-emission data base. 
Standard imaging and CLEAN (Hogbom 1974) proce- 
dures were then applied, and the resultant image was fur- 
ther improved by a standard self-calibration routine (e.g. 
Cornwell & Fomalont 1989). 

In the u-v domain, the continuum emission was re- 
moved from channels containing the H 1 data using 
the line- free channels (van Langevelde & Cotton 1990; 
Cornwell et al. 1992). After concatenation of the four 
Hi data sets, the gain solutions obtained from the self- 
calibration of the continuum image were applied. 

To facilitate subsequent analysis two sets of images 
were produced for both the continuum and line data. 
'Natural weighting' provided a resolution of 19" (R.A.) 
x 25" (Dec.) for investigation of the extended structure; 
'uniform weighting' yielded a restoring beam of 3.2" x 4" 
to reveal the spatial fine structure of the nuclear region. 
A correction for the gain variation across the ATCA 22-m 
antenna beams was applied to all images. 
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The rms noise of the final continuum image is ~ 
2.0 mJy beam -1 . Faint artifacts are present with flux den- 
sities of 5 - 8 mJy beam -1 . The peak-to-noise ratio for the 
nucleus is of the order of 800. Averaging two contiguous 
channels (this yields a channel spacing of 6.6 km s -1 ) the 
rms noise of the Hi data becomes 1.2 mJy beam -1 . 

To investigate the possibility of missing extended emis- 
sion, the galaxy was also mapped with the Parkes 64-m 
telescope. The total Hi flux density integrated over veloc- 
ity for Parkes observations (64 Jykms -1 ) was found to be 
9% lower than its ATCA counterpart (70 Jykms" 1 ). The 
difference is within the uncertainty of at least 5% for each 
estimate. In summary, the comparison of flux densities 
provides no evidence for missing emission in the ATCA 
observations. 

2.2. 12 COJ=2^1 data 

CO(2— 1) line emission was mapped during 1993 March 
1-10 with the SEST at an angular resolution of 23". A 
230 GHz SIS receiver yielded system temperatures of 600 
- 1300 K on an antenna temperature (T£) scale. The main 
beam efficiency was 0.46. An acousto-optical spectrome- 
ter, with 1440 channels and a total bandwidth of 1 GHz, 
provided a channel separation of 0.9 km s -1 . 

A rectangular grid of positions was selected, centred on 
the nucleus of the galaxy (see Table [I]), and with offsets 
parallel or orthogonal to the major axis at PA = 43°. The 
spectra were obtained at intervals of 10" along the major 
and minor axes, and 14" elsewhere. In total, 404 positions 
were sampled with offsets from the nucleus that ranged 
from —370" to +360" parallel to the major axis, and from 
—60" to +60" perpendicular to it. 

All spectra were obtained using a dual beam-switching 
mode (switching frequency 6 Hz), with a beam throw of 
11 '.7. The integration time was 4min per position and, 
averaging seven contiguous channels (channel spacing: 
6.3 km s -1 ), the rms noise ranged between 20 and 45 mK 
on a T£ scale. A set of three spectra was preceded by 
a short calibration observation of a black body 'paddle' 
which provided an intensity conversion to T£. Periodic 
continuum observations of the small-diameter nucleus of 
the radio source CenA at 115 GHz provided antenna 
pointing corrections. The pointing was also assessed by 
regularly re-observing the profile shape of the CO spec- 
trum at the central position of NGC4945. These mea- 
surements infer absolute positional uncertainties <8" and 
relative positional errors <5". 

The data were processed with the CLASS package of 
the Groupe dAstrophysique de Grenoble. To facilitate a 
comparison of CO with Hi, the CO spectra were trans- 
formed into a data cube and transferred to the AIPS soft- 
ware. 



3. Results and Discussion 




too fioo 800 




m 600 BOO 
VELOCITY [KM S" 1 ] 

Fig. 2. Line emission profiles as a function of Local 
Standard of Rest (LSR) velocity of (a) the Hi and (b) 
the CO(2— 1) emission integrated over the observed region 
with channel spacings of 6.6 and 6.3kms _1 , respectively. 

3.1. A brief overview 

The distributions of 1.4 GHz continuum emission, inte- 
grated Hi emission, and integrated CO(2— 1) emission are 
shown in Fig.|l|, superimposed on an optical image from 
the UK Schmidt SRC (Science Research Council) survey. 
The position offsets are relative to a(2000) = 13 h 05 m 
2754, £(2000) = -49° 28' 05". 

3.1.1. The radio continuum 

The continuum emission (Fig.[l]a) shows a bright small- 
diameter nucleus superimposed on an elongated distri- 
bution of emission. The map is in reasonable agreement 
with the 1.4 GHz image of Elmouttie et al. (1997) but 
shows additional details. The peak flux density of the 
central source is 4.2±0.1Jy beam -1 at a(2000) = 13 h 05 m 
27?3, 5(2000) = -49°28'07" (the position agrees to within 
2" with those given by Whiteoak & Bunton (1985) and 
Elmouttie et al. (1997)). The flux density integrated over 
the central source is 4.6+0. Uy. 
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400 600 KM S" 

Fig. 3. (a) H I spectrum (LSR velocities) towards the 
central continuum source, obtained with an angular reso- 
lution of 19" (R.A.) x 25" (Dec.) and a channel spacing 
of 6.6kms -1 . (b) The 12 CO(2— 1) line spectrum towards 
the centre position with an angular resolution of 23". The 
integration time was 60 minutes and the channel spacing 
is 0.9kms -1 . 



The emission extends over 11.6' x 3.3' with the ma- 
jor axis at a position angle PA ~ 45°. Northeast of the 
nucleus a ridge-line runs above the major axis and veers 
to the east at R ~ 240". Southwest of the nucleus the 
ridge appears to run below the major axis, turning to the 
west at R ~ 300". Faint (5 mJy) curved structures extend- 
ing out of the disk at R ~ 150-200" are residual sidelobes 
caused by the relatively bright emission at the nucleus (cf. 
Sect. 2.1). 




Offsets ["] 

Fig. 4. (a) Integrated CO(2— 1) emission convolved to 
43". Contours are 5, 10, 20, 30, 40, 50, 80, 160, 
320, 640 K beam" 1 on a T mh scale, (b) The CO 
(2— l)/CO(l— 0) ratio (in units of T m b) across the disk. 
The inner box denotes the approximate boundary for 'rea- 
sonably' convolved CO(2— 1) spectra. 

absorption was observed against the extended continuum 
outside this central region. Because of the influence of the 
absorption, the total integrated flux density discussed in 
Sect. 2.1 is not representative of the total Hi content. In 
the figure, the peak value of the integrated emission is 
~ 3.8 Jykms -1 beam" 1 (at R ~ +220"). The H I emission 
extends over ~ 22' x 4' (43 x 7.8 kpc), showing a stretched 
'S '-structure. Within this, a ridge of emission extends 13' 
across the nucleus at PA ~ 43°, whereas the major axis 
as defined by the 3% intensity level is at PA ~ 50°. At the 
south-western end of the distribution there is a moderately 
bright H I concentration that is extending to the north. At 
the opposite end of the galaxy the H I distribution extends 
to the south. 

Fig.|2pL shows the average Hi spectrum for the entire 
galaxy. The velocity relative to the Local Standard of Rest 
(LSR) extends from 340 to 770 kms -1 and the emission 
peaks at 400 and 715 kms -1 . The systemic velocity, as 
given by the average of the extreme and peak velocities, is 
V sys ~ 555 and 555-560 kms -1 , respectively (cf. Table |I]). 
The deep asymmetric central depression reflects the pres- 
ence of absorption against the nucleus. The shape of the 
spectrum is similar to that obtained with the Parkes tele- 
scope directed towards the centre of the galaxy (Whiteoak 
& Gardner 1977). 

Fig.|3^ shows our Hi spectrum with ~23" resolution 
towards the position of the nuclear continuum source. Hi 
is present in absorption against the continuum emission at 
velocities of 350 - 770 kms -1 , i.e. symmetrically offset up 
to ^210 kms -1 from V sys . The profile is asymmetric and 
appears to be composed of several components. A peak 
absorption flux density of 1.05 Jy beam -1 is observed at 
570 and 635 kms -1 . The corresponding linc-to-continuum 
ratio is 0.25. 



3.1.2. The Hi distribution 

In our contour map of the Hi distribution (Fig.|l|b) the 
central region is dominated by H I absorption against the 
strong continuum emission of the nucleus. However, no H I 



3.1.3. CO emission 

The distribution of CO(2— 1) emission integrated over ve- 
locity (Fig.|l|c) extends over more than 10' along the ma- 
jor axis. As the continuum and H I distributions, CO also 
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shows a ridge extending above the major axis north-east 
of the nucleus and below the major axis to the south-west. 
Like the continuum the CO distribution bends towards the 
east at R ~ 250". 

Fig. ^|b shows the integrated CO(2— 1) spectrum for the 
entire region observed. The CO velocity range is about the 
same as for the H I in Fig.|^a. Broad peaks are present near 
415, 480, 580, and 705 km s -1 . For an assumed symmetric 
system, the outer peaks would suggest a systemic velocity 
of 560 km s -1 . 

Fig.||b shows the CO(2— 1) emission towards the nu- 
cleus of the galaxy. The CO covers the same velocity range 
as the previous spectrum. The profile shape is similar, ex- 
cept that the higher velocity feature is fainter than the 
systemic feature. In our Fig.^j and in Fig. 5 of Dahlem et 
al. (1993), both features have almost the same line temper- 
ature. This likely reflects small differences (a few arcsec) 
in the pointing of the telescope (cf. Sect. 2.2). A fitting of 
gaussian components yields distinct components centred 
at velocities of 447, 493, 593, and 701 km s -1 , with a fur- 
ther underlying broad component centred at ~565kms _ . 

Fig.^a shows the integrated CO(2— 1) emission con- 
volved to a resolution of 43". Combined with corre- 
sponding SEST CO(l— 0) data observed at the same 
resolution (Dahlem et al. 1993), the distribution of the 
CO(2-1)/CO(1-0) ratio is shown in Fig.@b. The ratio 
varies from 0.8 to 2.0 and demonstrates that 'warm spots' 
with ratios larger than unity are not confined to the cen- 
tral region but are also observed far out in the disk. For 
a possible spatial correlation of these warm molecular re- 
gions with spiral arms, see Sect. 3.4. 

3.2. The central region 

The strong Hi absorption against the central radio con- 
tinuum source complicates a direct comparison of H I and 
CO. The Hi absorption must originate from in front of 
the continuum whereas the CO emission may arise from 
in front of and behind the nucleus. This difference is con- 
sistent with the CO and Hi lineshapes shown in Fig.||, 
where the CO profile is wider than its Hi counterpart at 
the half-maximum-intensity points. This effect must be 
significant: We could also plot, instead of the Hi absorb- 
ing flux, the Hi optical depth r(Hi) = -ln(l + S L /FS C ) 
with F denoting the continuum source covering factor and 
5l and S c being the (negative) line and (positive) contin- 
uum flux. For F=l the profile would resemble that shown 
in Fig. ||a, since the line remains optically thin even at the 
line centre. With F ~ 0.25, however, Hi optical depths 
would be large near V ~ 600 km s" 1 and the Hi column 
density profile would become narrower. The velocity of 
the central CO peak lies between the velocities of the two 
strongest H I absorption components. 

Fig. ||a shows the nuclear continuum emission at a res- 
olution of 3 '.'2 in R.A. and 4"0 in Dec. It is spatially 
resolved: a gaussian fit and a correction for beam size 
yields a source size of 7 '.'6 x 3 '.'4 (±0"2) (correspond- 




: 1 15 5 • -5" -15" ' 
OFFSET ALONG MAJOR AXIS 

Fig. 5. The resolved nucleus at a resolution of 3.2" in 
R.A. and 4'.'0 in Dec. (a) The nuclear 1.4 GHz contin- 
uum. The contour levels have flux densities of -2 (dashed) , 
10, 20, 40, 60, 120, 240, 500, 1000, 1300 mJy beam- 1 , 
(b) The LSR velocity field of the nuclear Hi absorp- 
tion with iso-velocity contours of 435, 460, 485, ... 685, 
and 710 km s -1 . The vertical direction is equivalent to 
PA = 43°. (c) Position- velocity diagram of the nuclear 
Hi absorption along PA = 43° at ~3.6" x 7.25 km s -1 
resolution. The velocity axis is labeled with respect to 
560 kms" 1 . Contours are -1, 1, 2,5, 10, 20, 50, 70, 85, 90, 
95, 99% of the peak absorption of 612 mJy (for the peak 
flux density in a ~23" beam, see Fig.^a). 
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ing to 250 x llOpc), with the major axis at PA = 42?5 
and a peak intensity of 1.33±0.03mJy beam -1 . The angu- 
lar dimensions are reasonably consistent with those (5 '.'7 
x 2"0 (±0'.'l)) determined by Whiteoak & Wilson (1990) 
at 6 GHz. Figs. ^> and c show the velocity field of the 
H I absorption towards the resolved nucleus and the as- 
sociated position-velocity diagram. The velocity field is 
centred within 15kms _1 of the galaxy systemic velocity 
(555-560 km s _1 ) and the velocity gradient across the nu- 
cleus is reasonably uniform. It supports the existence of 
a cloud which surrounds the nucleus and has rapid solid- 
body rotation. These results are consistent with the 6 GHz 
OH study by Whiteoak & Wilson (1990). However, be- 
cause the H I absorption is associated only with the region 
overlying the continuum emission, it cannot provide infor- 
mation about the overall size of the cloud. 

The strong CO(2— 1) emission associated with the nu- 
cleus (cf. Fig.|l|c) has an elongated distribution. Its extent, 
derived from a gaussian fit and corrected for beam size, is 
37" x 21" with PA = 13°. Table| contains a comparison 
with results for CO(l— 0) emission (Dahlem et al. 1993) 
and CO (3— 2) emission (Mauersberger et al. 1996a). The 
emission associated with the highest-state J = 3-2 tran- 
sition is more centrally concentrated than for the lower- 
excitation transitions. 

Fig.|| shows the CO(2— 1) position-velocity diagram 
taken along the major axis. Like Fig.|^c, it supports the 
presence of a molecular cloud centred on the nucleus and 
rotating as a solid body with a rotational velocity that 
reaches V ~ 160 km s . The difference between | V — ^^sys | 
~ 160 km s -1 and the edge of the Hi absorption, | — V^ ys | 
~ 210kms~ 1 (Sect. 3.1.2), is likely reflecting the inter- 
nal lincwidth of molecular clouds, the presence of excen- 
tric orbits, radial motions caused by cloud collisions, or 
non-axisymmetric distortions of the gravitational poten- 
tial (see e.g. Mauersberger et al. 1996b and Sect. 3.5.5). 
Radius (7 "5; taken from Fig.||) and rotational velocity 
(~160kms _1 ) are compatible with those obtained for 
the CO(l— 0) and CO(3— 2) rings observed by Dahlem et 
al. (1993) and Mauersberger et al. (1996a), respectively. 
Bergman et al. (1992) were the first to model the nuclear 
ring. Henkel et al. (1994) and Curran et al. (2001) com- 
mented that nuclear gas densities significantly higher than 
the originally proposed value (150cm -3 ) are required. 
Furthermore, the central CO(2— l)/CO(l— 0) ratio is a 
factor of two higher than that presumed by Bergman et 
al. (1992). 

The position- velocity diagrams (Figs.[|c and ||) show 
that for R < 6", the cloud rotates like a solid body, but at 
R > 10" the rotational velocity decreases progressively. If 
this velocity variation is interpreted in terms of rotation 
around a central mass confined essentially within R = 7 "5, 
then the estimated mass is 1.5 10 9 M Q . This is consistent 
with a value of 0.8-1.6 10 9 M Q within a radius of R ~ 
3", as suggested by Koornneef (1993). Mauersberger et 
al. (1996a) estimated a dynamical mass that, scaled to 
a galaxy distance of 6.7 Mpc, increases from 1.410 9 M Q 
within R = 5 "5 to -5 10 9 M n within i?=33". 



Table 3. Extent of the nuclear CO emission. A gaussian 
beam shape and a 2-dimensional gaussian approximation 
for the molecular cloud distribution have been assumed. 



12 CO 


Apparent 


Beam 


Deconvolved 


Ref. 


line 


FWHM 












["] 


["] 


["] 


[pc] 




(1-0) 


54±2 


43 


33±4 


1100±130 


1 


(2-1) 


32±2 


24 


21±4 


680±130 


1 




44±1 


23 


37±4 


1210±30 


2 (a) 




31±1 


23 


21±4 


680±30 


2 < b > 


(3-2) 


22±2 


15 


12±3 


360±130 


3 



(1) Dahlem et al. (1993), 

(2) this work; major and ^ minor axis values with respect 
to PA = 13°, 

(3) Mauersberger et al. (1996a), adjusted for D — 6.7 Mpc 



Fig.|7p.,b shows Hi position-velocity diagrams at two 
position angles (49° and 43°) to best represent the Hi 
motions along the major axis. A noteworthy feature is 
Hi emission from the central region at velocities above 
>200kms -1 relative to V sys: and higher than the veloci- 
ties seen in H I absorption or CO emission in Fig. | This 
Hi emission (there is no counterpart at V — V sys < - 
200 km s -1 ) may be related to an optical outflow from the 
centre region towards the halo, first suggested by Nakai 
(1989) on the basis of optical observations. Spectra of 
the Ha and [Nil] line emission close to the nucleus show 
high velocities consistent with a conical outflow of gas 
(Heckman et al. 1990). 

3.3. The overall gas distribution 

Qualitatively, it appears that, in spiral galaxies, neu- 
tral atomic gas dominates in the outer regions, whereas 
molecular gas is concentrated towards the inner disk (e.g. 
Sanders et al. 1984). Generally, the CO(l— 0) line integral 
is used to estimate the molecular mass content, conven- 
tionally expressed as the H2 mass. For NGC 4945, we have 
used the CO(2— 1) data on the basis that the uncertainty 
due to variations in the CO(2-1)/CO(1-0) ratio (Fig.|b) 
is considerably less than the uncertainty in the CO to H2 
conversion factor X. 

Velocity- integrated H 1 and CO line intensities were de- 
rived from spectra taken along the major axis between 
r = 100" and 360". The data indicate that the Hi column 
density increases, and the CO column density decreases, 
with radius. Using the 'standard' conversion factor X — 
2 x 10 20 cm _2 /(Kkms _1 ) we find for the inner region, at 
e.g. r ~ 135", that the column density ratio of H1/H2 is 
of order 0.2. At the edge of the observed CO(2— 1) disk, 
i.e. at R ~ 320 - 350", the H1/H2 column density ratio 
becomes ~1. However, these ratios should be treated with 
caution, because Mauersberger et al. (1996a, b) presented 
compelling evidence that, at least for the central regions 
of NGC 253 and NGC 4945, the value of X is almost an 
order of magnitude lower than the standard value. Thus 
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OFFSET ALONG hAJOR AXIS 

Fig. 6. CO(2-l) position-velocity diagram of NGC4945 
along the major axis (PA = 43°) with a resolution of 
23" and a channel spacing of 6.3 km s -1 . The velocity axis 
is labeled with respect to 560 km s -1 . Contours are 75, 
150, 300, 600, 1200, 1600, 1800, 1900 mK on a T mb scale. 
The rms noise is ^50 mK. For the corresponding CO(l— 0) 
diagram, see Fig. 4 of Dahlem et al. (1993). 

the H 1/H2 abundance ratio should also be ~1 in the inner 
region. 

The velocity-integrated H 1 emission (Fig. ^Jd) suggests 
a variation in major axis position angle and perhaps incli- 
nation across the galaxy. As mentioned earlier, there ap- 
pears to be an H 1 inner disk with a lower position angle 
than that for the outer H 1 regions. To investigate whether 
the galaxy is warped, we modeled the H 1 emission with 
a set of concentric rings (cf. Begeman 1987; Rogstad et 
al. 1974) of width 50" from R = 100" to 650" and fitted 
the north-eastern and south-western sides of the galaxy 
separately. The resulting position angles and inclinations 
for the rings both vary over less than 10° along the ma- 
jor axis. Part of this variation is systematic in that the 
position angle varies from 50° in the north-east to 43° in 
the south-west, and for |i?|>400" the inclination appears 
about 3° higher in the south-west than in the north-east. 
Significant variations in inclination are also found as a 
function of azimuth and radius in M31 (see Braun 1991) 
and may be in part attributable to the presence of spiral 
arms. 

3.4. Spiral arms 

It has been pointed out already that the distributions of 
the 1.4 GHz continuum, Hi, and CO in Fig.[l], and also 
the optical image (for optical images not contaminated by 



radio contours, see e.g. Nakai 1989; Elmouttie et al. 1997) 
show a common extended emission ridge which is above 
the major axis at positive offsets, and below the major axis 
at negative offsets. There are many emission features as- 
sociated with the ridge, and the most prominent are listed 
in Table ||. Several emission peaks are located at common 
major axis offsets, and are presumably related physically. 
The H 1 distribution shows distinctive features at each end 
of the galaxy, at offsets of +620" and —580". They extend 
away from the major axis and resemble the tangential lo- 
cations of outer trailing spiral arms of a highly inclined 
galaxy. At the north-eastern end of the galaxy, which is 
moving away from us, the arm ending at the near side ex- 
tends around behind the plane, above the major axis, and 
joins an inner disk at the location of the emission features 
at offset ~ —250". At the south-western end, rotating to- 
wards us, the arm extends around the front of the galaxy, 
below the major axis, joining the inner disk at the loca- 
tion of the emission peak at offset +230". The prominent 
dust lane present in the optical image would be associated 
with this arm. This interpretation accounts for the anti- 
symmetry of the inner emission ridge (see also Dahlem et 
al. 1993 for the inner part of the galaxy). 

It is noteworthy that strong IRAS point sources are 
located at R ~ 240" on each side of the nucleus. At these 
locations the 12 CO(2-1)/ 12 CO(1-0) line intensity ratio is 
>1. These features are consistent with regions of enhanced 
star formation, which may reside in spiral arms viewed 
near their tangential point. 

The interpretation could be extended to include the 
regions of enhanced emission at offsets of about -120" 
and +120" as additional tangential locations of the two 
spiral arms. Then one arm would extend through offsets 
+620", -250" and +120", and the other through offsets 
-580", +230" and -120". 

Interpreting the distributions in Fig. [l] in terms of a 
two-arm spiral structure, we can test whether the desig- 
nated tangential locations are consistent with an overall 
regular spiral pattern. We therefore consider a spiral de- 
scribed by R — Rq x exp a?9, with the pitch angle tp in 
degrees given by a = tan tp and the azimuth angle $ in ra- 
dians. Consideration of the two inner tangential offsets for 
each arm yields tp — 12° and 13°. For a pair of positions 
consisting of the outermost and central offsets, tp = 16° 
for both arms. The results are plausible - two major spiral 
arms with pitch angles of about 12 — 13° in the inner re- 
gion and 16° in the outer region. They are consistent with 
arms in other spiral galaxies (e.g. Puerari & Dottori 1992) 
and might also be detectable by near infrared imaging. 



3.5. Disk kinematics and mass 
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Table 4. Positions of regions with enhanced continuum, 
H I or CO emission. Positional errors are <25". 

Offsets in arcsec 



Cont. Hi CO(2-l) 





+620 




+230 


+230 


+230 






+170 


+130 


+120 


+90 


-55 




-80 


-110 


-120 


-140 


-175 


-185 


-195 


-215 


-235 


-250 


-255 


-270 




-290 


-320 






-335 


-350 




-580 





3.5.1. Modelling the overall velocity distribution 

A commonly used relation to describe the rotation curve 
of a galaxy is 



Vmax (| + |x(^-)«)(3/2«) 

(see Eqs. 26 and 28 of Brandt et al. 1960), where n is a 
measure of the steepness of the rotation curve and i?v max 
is the radius at which the maximum rotation velocity V max 
occurs. This curve results in solid-body rotation at small 
radii and a Keplerian-type velocity decrease at large radii. 
For a flat spiral galaxy inclined to the line-of-sight (such 
as NGC4945), the rotation curve can be derived from the 
variation of radial velocity along the major axis, assuming 
that the maximum velocity at any position, corrected for 
the galaxy inclination, represents the tangential (i.e. cir- 
cular) velocity at the corresponding galactocentric radius. 

Inspection of the H I position- velocity diagrams along 
the major axis (Figs. 0a, b) suggests the existence of two 
general velocity regimes. The rapidly rotating central 
molecular cloud has been discussed already in Sect. 3.2. 
The second system relates to the main disk of the galaxy, 
viewed at a position angle of ~ 45°. It extends out to R ~ 
400", with rotational velocities reaching V ~ 170 km s -1 
(the correction for an inclination of 78° is only 2% and is 
ignored) . This disk is associated with the inner concentra- 
tion seen in the distribution of integrated Hi (Fig.[l]b). In 
the position- velocity diagrams the additional H I concen- 
trations at R ~ 600" then correspond to the spiral arms 
that have already been discussed (Sect. 3.4). 

Because of the high inclination of the galaxy, a line- 
of-sight at a specific offset along the major axis will also 
include H I at galactocentric radii larger than the radius 
corresponding to that offset. The line-of-sight velocities of 
the extra H I components are smaller than the rotational 
velocity at a specific offset and are responsible for the low- 
velocity 'tail' present at many positions along the major 
axis. 



Table 5. Parameters derived from the Brandt-curve fits 
and used to model the large scale velocity field for the H I 
disk. 



Hi 



Kinematical centre 


a(2000) 


13 h 05 m 25f34 




5(2000) 


-49°29'09" 


Systemic velocity 




557±3kms- 1 


Max. rotation vel. 


Vmax 


170±5kms- 1 


Turn-over point 


-^Vmax 


380" ± 10 


Position angle 


PA 


45° ± 2 


Inclination angle 


i 


78° ± 1 


Brandt curve index 


n 


3.5+1 



Single gaussian profiles were fitted to individual spec- 
tra in the Hi dataset. To avoid the effects of wings in 
the H I spectra, only the top half of the H I profiles with 
Sh > 5'peak/2 was used in the fitting. Applying an iter- 
ative procedure, Eq. 1 was fitted to the profile velocities 
for offsets out to ±450". The parameters listed in Table 
H were obtained for the disk region. They are consistent 
with the values obtained from a similar fit to CO(l— 0) 
data by Dahlem et al. (1993). Our n values are 3 (for Hi) 
and, from a similar fit to the CO(2— 1) data, 4, with small 
formal errors. Our value for n is significantly smaller than 
that determined by Abies et al. (1987), n ~ 7. This may 
be a consequence of the better angular resolution of our 
data. Inclusion of the outermost features in the analysis 
would have increased -Rv max ; this value is not well-defined 
in the fit to the H I position- velocity diagram and the small 
uncertainty in Table ^ only accounts for the formal errror. 

3.5.2. The mass of NGC4945 

The mass of NGC 4945 can be derived with 

M = (1.5) * x 1242 xl^xJUx D (2) 

from the parameters in Table |E| (see e.g. Eq. 29 of Brandt 
1960; Rogstad et al. 1967; Dahlem et al. 1993). M is in 
solar masses, V max is in kins" 1 , i? m ax is in arcsec, and D 
is in Mpc. For V max = 170 km s" 1 and i? max = 380", M 
= 1.4 10 11 Mm. Reference to the Hi position-velocity di- 
agrams (Fig. [7]) suggests that this value may be a lower 
limit, because the velocities of the outermost Hi features 
are larger than the values given by the fitted rotation 
curve. The paucity of Hi at offsets of 400 — 500" may 
signify a lack of gas at the tangential location along the 
lincs-of-sight, in which case the observed velocities are less 
than the rotation velocities (see also Sect. 3.5.6). 

The total integrated Hi flux density of 70Jykms~ 1 
yields an Hi mass of 7.5 x 10 s M Q , which is only 0.5% 
of the total mass and considerably less than expected for 
such galaxies. However, the presence of the central H I ab- 
sorption has badly affected the integration. Exclusion of 
the H I spectrum observed towards the peak of the nuclear 
continuum emission (see Fig. ^a) would have increased the 
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integrated flux density by more than 200 Jy km s , yield- 
ing a total mass fraction of 2%. 

Because of the limited coverage of the CO distribu- 
tion, the molecular content of the galaxy was not esti- 
mated from our CO(2— 1) data. Assuming an integrated 
CO intensity to H2 column density conversion factor of 
X = N(R 2 )/Ico = 2.3 1 20 cm- 2 (Kkms -1 )- 1 , Dahlem 
et al. (1993) derived a molecular gas mass fraction of 5% 
from their CO(l— 0) data extending out to 360". While 
the X-factor is appropriate for the local Galactic disk, it 
may be too large by a factor of 2-5 for a galaxy with 
a strong nuclear contribution to the overall CO emission 
(see e.g. Mauersberger et al. 1996a). Thus Hi and molec- 
ular masses in NGC4945 should be similar, at the order 
of 2% of the total mass inside a radius of 380" (12 kpc). 

3.5.3. Kinematical fine structure 

To compare the H 1 and CO kinematics on a small scale, 
pairs of spectra were taken at offsets between —100" and 
—360" along the major axis, but with a 10" offset to 
the south; these directions provided the highest CO in- 
tensities for comparison. Velocity differences between Hi 
and CO(2— 1) features were found to range from 5 to 
20kms~ . They are not systematic and show a random 
variation in sign. Aside from the existence of H 1 and CO 
clouds with separate peculiar motions, for a highly in- 
clined galaxy such as NGC 4945, differences may arise be- 
cause a number of clouds at different physical locations 
within the galaxy are seen simultaneously in one beam. 
A line-of-sight in a specific direction along the major axis 
may include clouds located at different galactocentric radii 
or even at different heights above and below the galactic 
plane. 

In the first case (different galactocentric radii), a line- 
of-sight displacement of 1.5 to 3 kpc between Hi and CO 
could explain velocity differences of 5 to 20kms~ 1 . For 
comparison, the CO emission peaks in M51 are typically 
separated by 1.5 to 2 kpc along the spiral arms (Garcia- 
Burillo et al. 1993). Much smaller line-of-sight displace- 
ments could explain the observed velocity differences if 
streaming motions due to density waves were prominent 
(for M51, see Aalto et al. 1999). 

In the second case (different heights above the plane), 
significantly slower rotation velocities might exist in the 
halo compared to the underlying disk, as observed in the 
case of M82 (Sofue et al. 1992). Also in our Galaxy pe- 
culiar velocities are observed in high-latitude H 1 gas (e.g. 
de Boer 1985; their Fig. 2). 

3.5.4. Velocity fields and velocity residuals 

Fig. ||a shows the Hi velocity field. The distribution shows 
anti-symmetric contour distortions with respect to the ma- 
jor axis, indicating departures from uniform circular mo- 
tion. Considering a given velocity isophote north-west of 
the major axis, contours are displaced towards the north, 



whereas south-east of the major axis contours are dis- 
placed to the south. The corresponding CO (2-1) velocity 
field indicates a similarly anti-symmetric behaviour. 

Fig. [pa displays the model velocity field obtained using 
Eq. 1 and the parameters listed in Table |5|. To examine the 
departure of H 1 motions from the estimated circular ro- 
tation, the model field was subtracted from the observed 
Hi field. Fig.||c shows the distribution of Hi residuals. 
In the outer parts of the galaxy the residuals are typ- 
ically 0-20 km s -1 with no systematic trends, indicating 
that motions are to first order consistent with uniform 
circular motion. However, closer to the nucleus, i.e. above 
the major axis of the galaxy in the north-east and below 
the major axis in the south-west, higher residuals of up to 
±40kms~ 1 define an elongated region extending across 
the nucleus to offsets of +120" and -200". The system- 
atic departures in H 1 velocities from circular motion were 
first identified by Abies et al. (1987), but the lower resolu- 
tion of their Hi data (47"x73") led to a more widespread 
contamination by the nuclear absorption and inhibited the 
detection of velocity departures as large as 40kms~ 1 . 

Abies et al. (1987) detected velocity residuals of 15 to 
30kms _1 south-west and east and of -15 to -30kms _1 
west and north-east of the nucleus. This was interpreted 
in terms of radial motion towards the nucleus. The resid- 
uals of our fit are slightly different both with respect to 
amplitude (as already mentioned) and geometry: In the 
north-east, north of the major axis, we only detect signif- 
icant negative velocity residuals, while in the south-west, 
south of the major axis, large positive residuals are found 
(Figjc). 

Our Fig.|| shows a prominent Hi feature displaced by 
about +80kms _1 from the systemic velocity. Since this 
feature is seen in absorption against the nuclear contin- 
uum source, it signifies neutral atomic gas approaching 
the galactic centre. From their CO(l— 0) and (2—1) data, 
Dahlem et al. (1993) find evidence for inflowing low den- 
sity molecular gas at velocities of 80kms _1 . It seems that 
our nuclear Hi absorption spectrum that separates the 
systemic from the inflowing velocity component (for a less 
conclusive spectrum, see Abies et al. 1987) is tracing the 
neutral atomic part of this kinematical feature as far as it 
is located in front of the nucleus. Note that this gas com- 
ponent is located inside the nuclear molecular ring (see 
Table |; Bergman et al. 1992; Dahlem et al. 1993), at a 
galactocentric radius of a few hundred pc or less. 

3.5.5. Evidence for a bar 

De Vaucouleurs (1964) suggested that NGC 4945 contains 
a central bar, which provides a very efficient mechanism to 
quickly transport matter towards the nucleus while releas- 
ing angular momentum outwards. Can such a bar explain 
the velocity anomalies outlined in Sect. 3.5.4? A classical 
signature of a bar is an S-shaped distortion in the isove- 
locity contours of the gas (e.g. Kalnajs 1978; Duval & 
Monnet 1985). This is observed in NGC4945 (see Fig.|a; 
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Sect. 3.5.4). If a bar is present, the Hi velocity residu- 
als in Fig. ||c could be interpreted in terms of systematic 
overall gas flows along it, approaching at positive offsets, 
receding at negative offsets. The orbit velocities in the 
bar near the nucleus may be estimated from the pv dia- 
grams and the velocity field residuals. The continuum and 
CO distributions in Figs.|l| and |] show enhanced emis- 
sion extending out from the nucleus along PA ~ 33°, i.e. 
along the axis that is also showing the velocity anoma- 
lies (Fig.||c). Assuming that the putative bar is associated 
with the emission ridge introduced above and that this is 
located within the plane of the galaxy, we obtain an az- 
imuthal angle of ~ 40° (counterclockwise) with respect to 
the line-of-sight projected onto the plane of the galaxy. 
The total inclination to the line-of-sight would then be 
45°. 

There exists a strong correlation between the nuclear 
absorbing column density and the presence of a bar in 
Seyfert 2 galaxies. Strongly barred Seyfert 2 galaxies have 
an average iVjj that is two orders of magnitude higher than 
in non-barred Sy 2s. More than 80% of the 'Compton thick 
Seyfert 2s' (JVh <; 10 24 cm -2 ; most of this column density 
must arise from the innermost few 10 pc) are barred and 
almost 60% of these are 'strongly' barred (Maiolino et al. 
1999). NGC 4945 is a Compton thick Seyfert 2 galaxy (e.g. 
Guainazzi et al. 2000; Madejski et al. 2000). 

Non-interacting spiral galaxies with moderate inclina- 
tion and substantial far infrared emission (for the cen- 
tral part of NGC 4945, L F IR~ 5 10 10 L ; IRAS 1989) are 
known to show both strong and long bars (Martinet & 
Friedli 1997). The presence of asymmetric morphologies 
in individual Seyfert galaxies is positively correlated with 
their tendency to exhibit enhanced star forming activity 
(Maiolino et al. 1997). Hence irrespective of direct obser- 
vational evidence, the presence of a bar in NGC 4945 is, 
at least statistically, very likely. 

The velocity dispersion of H I in the central region is 
also consistent with the presence of a bar. Estimated for 
each observed Hi spectrum during the gaussian fitting, 
the disperson in the outer galaxy is 10 - 30kms _1 , which 
is consistent with the dispersion of gas along the spiral 
arms in the presence of relatively uniform circular rota- 
tion. However, within ~ 150" of the nucleus there is an 
elongated region where the dispersion is above 50kms _1 . 
This may reflect higher gas turbulence or a higher space 
density of gas clouds near the centre, but could also reflect 
fast motion on highly eccentric orbits in a bar. Since all 
this evidence is circumstantial, a definite answer to the 
question whether NGC 4945 has a bar must come from 
the near-infrared (preferably K-band) image. This method 
proved to be conclusive in the case of another nearby 
highly inclined southern starburst galaxy, NGC 253 (see 
e.g. Engelbracht et al. 1998). 



3.5.6. Dynamical analysis 

The Hi and CO(2— 1) pv diagrams can be analysed using 
linear resonance theory (e.g. Binney & Tremaine 1987) to 
deduce the locations of various gravitational resonances 
within NGC 4945. If we assume the presence of a 'weak' 
bar (the orbits can be described by a superposition of 
circular motion around the centre and small oscillations 
caused by the non-axisymmetric part of the potential), 
the spiral and barred structure is constrained by the loca- 
tions of these resonances. For a perfectly edge-on barred 
spiral the locations can be only approximate, because the 
rotation curve will be affected by the resonances as well as 
the bar's de-projected size and angle to the line-of-sight. 
In the case of NGC 4945, the proposed bar (Sect. 3.5.5) 
is sufficiently displaced from the major axis to provide a 
reasonable approximation. 

For gas orbiting the centre of a galaxy with angular ve- 
locity Q = V/R at radius R, the radial epicyclic frequency 
k is expressed by 



K (R) = \ Un 2 (R) + R 



dn 2 (R) 

dR 



(3) 



Application of a Brandt rotation curve (Eq. 1) infers 

n-l 



dfl 2 
~dR 



2 V 2 



R 



i?Vma 3 



2 

3 \i?v 



R 



(4) 



Inner and outer 'Lindblad' resonances (ILR and OLR) oc- 
cur when the pattern speed of the bar Qp = Q — n/m and 
f2 + k/tti, respectively (m—1 represents a barred poten- 
tial, see e.g. Eq. 3-115 of Binney & Tremaine 1987). Fig.| 
shows the measured variation of angular velocity as a func- 
tion of galactocentric radius. A rotation curve of the form 
given in Eq. 1 with the parameters outlined in Table || has 
been fitted to Hi velocities derived from the Hi pv dia- 
gram m Fig.g. The fitted curve provides a reasonable fit 
to the outer galaxy (R <; 5 kpc) . 

For a barred spiral galaxy, it is expected that the 
co-rotation radius CR would be located not far beyond 
the end of the bar. For NGC 4945 this position is un- 
certain. Accounting for unsystematic velocity residuals of 
- 20kms~ 1 (Sect. 3.5.4), the velocity anomaly can be 
traced out to approximately ±150" projected on the sky 
(see Fig.||c). With an inclination to the line-of-sight of 
45° (Sect. 3.5.5) this would be equivalent to a bar extent 
(from the nucleus) of about 7 kpc. Since the spiral arms 
could be traced from the outer galaxy (^600" offset from 
the centre) to ~120" (Sect. 3.4), the spiral arms might be 
connected with the bar, as commonly observed in galax- 
ies with small inclination (e.g. Reynaud & Downes 1997; 
Hiittemeister et al. 1999). The CR must be located at 
f?>7kpc, the outer radius of the gaseous bar. Since stellar 
bars are longer than their gaseous counterparts (Martinet 
1995) and since gravitational torques de-populate the co- 
rotation region in a spiral galaxy (e.g. Garcfa-Burillo & 
Guelin 1995; Combes 1996), the drop in Hi intensity seen 
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Fig. 9. Variation of angular velocity (f2) as a function of 
galactocentric radius, determined from Hi observations 
(asterisks; see Sect. 3.5.1 for details) and fitted with a 
Brandt rotation curve (solid line; see Eq. 1) for the outer 
disk with parameters given in Table §. At R ~ 15 kpc, 
low measured angular velocities with respect to the fit are 
likely an artefact caused by a lack of neutral atomic gas at 
this radius. Shown are also fl—n/m and Q+n/m for m = 2 
(dashed lines) and m — 4 (dotted lines). The assumed co- 
rotation radius (vertical dashed line) and resulting pattern 
speed rip (horizontal dashed line) are indicated. 

in Fig. [7] at R ~ 450" (14.6 kpc) and the small 'rota- 
tional' velocities at this radius (see Fig. ^| and Sect. 3.5.2) 
may be considered as a signature of the CR (see also 
Freeman 1997 for 'average' radii). Adopting this radius, 
the pattern speed is llkms -1 kpc -1 ; the OLR (which oc- 
curs where the pattern speed intercepts the upper m=2 
curve in Fig.^) is then at a radius of 21 kpc (angular dis- 
tance: ~650") and the m=4 ultra-harmonic resonance is 
at 18 kpc (550"), where the outermost H 1 features are seen 
(Fig. U). 

Our Brandt rotation curve is a good fit to the outer 
galaxy, but does not allow to predict the presence of an 
ILR. In the innermost parts of the galaxy, measured angu- 
lar velocities are larger than those suggested by the Brandt 
curve (Fig.^), leaving open the possibility that an inner 
Lindblad resonance exists at R < 3 kpc. An attractive but 
speculative view is to associate the inner molecular ring 
at a galactocentric radius of a few hundred pc (Table ||) 
with the ILR that might contain a nested secondary bar 
(the 80kms _1 inflow; see Sect. 3.5.4) guiding atomic and 
molecular gas to the putative circumnuclear torus discov- 
ered by Greenhill et al. (1997). Note however that the 
nuclear molecular ring is part of the transitional region 



between solid body and differential rotation. The ring is 
therefore not necessarily formed by a gravitational res- 
onance. It could also be caused by viscous transport as 
e.g. outlined by Dather & Biermann (1990). Adopting the 
Dather & Biermann mechanism, an age estimate applying 
their last equation leads to a formation timescale of i r i ng 
-lGyr. 

4. Conclusions 

A high-resolution study has been performed of the 1.4 GHz 
continuum, Hi, and CO(2— 1) emission for the southern 
spiral galaxy NGC 4945. It utilizes both the Australia 
Telescope Compact Array (ATCA) and the Swedish-ESO- 
Submillimetre Telescope (SEST). The angular resolution 
is ~23" (750 pc at Z?=6.7Mpc) and the spectral resolution 
is ~7kms _1 . The ATCA results also yield high resolution 
(~3.6") images of the nuclear region of the galaxy. 
The main conclusions are as follows: 

1. The 1.4 GHz continuum emission extends over 
11.6'x3.3' with PA ~ 45° and contains a bright nu- 
cleus with a size of 7'.'6x3"4 (±0"2) (corresponding 
to 250xll0pc) with the major axis at PA = 42?5 and 
an integrated flux density of 4.6 ± O.lJy, centred at 
a(2000) = 13 h 05 m 2753, J(2000) = -49°28'07". 

2. Our CO(2— 1) data confirm the presence of a circum- 
nuclear ring at R ~ 10-20". H 1 gas is observed near the 
centre with a rotation velocity that reaches 160 km s -1 
at R=7'.'5 (250 pc); at smaller and larger radii the ve- 
locity is lower. The dynamical mass inside of 7 "5 is 
1.5 10 9 M Q . 

3. Nuclear Hi emission at velocities offset by more than 
200 km s -1 from the systemic velocity may be related 
to an outflow towards the halo, observed at optical 
wavelengths. Hi absorption shows evidence for inflow- 
ing nuclear gas (|V — V sys \ ~ +80kms _1 ) that was 
previously only seen in CO. 

4. H 1 is detected beyond the optical 25 m arcsec -2 con- 
tours. The Hi emission extends over ~ 22' (43 kpc) at 
a position angle that increases from 43° to 50° from 
the south-west to the north-east. Hi Local Standard of 
Rest velocities extend from 340 to 770 km s^ 1 and yield 
a systemic velocity of 555-560 km s" 1 . An analysis of 
the velocity field indicates that the maximum mean ro- 
tation velocity of 170 km s -1 occurs at a galactocentric 
radius of ^380", yielding a total mass of 1.410 n M Q 
inside of R — 380". Mass fractions of molecular and 
Hi gas are ~2% for each component. 

5. The Hi distribution shows distinctive features at each 
end of the major axis, offset about ±600" from the 
nucleus. These are interpreted as spiral arms viewed 
tangentially. Other common features in the continuum, 
Hi and CO distributions at offsets near ±120" and 
±240" may be additional tangential directions for the 
arms. The results are consistent with two spiral arms 
having pitch angles of 12-13° in the inner (|i?|<240") 
and 16° in the outer galaxy. 
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6. The H I velocity field of the galaxy shows a pattern gen- 
erally consistent with uniform circular motion. Near 
the nucleus, however, systematic departures exceeding 
±30kms _1 , S-shaped velocity contours, high infrared 
luminosities, large nuclear gas column densities, and 
increased velocity dispersions are interpreted in terms 
of a bar extending across the nucleus. The data sug- 
gest a structure with PA ~ 35° , an azimuthal angle of 
~40° (counterclockwise) , extending to radial offsets of 
±150" and connecting the inner spiral arms. The likely 
presence of a bar and large amounts of molecular gas 
in the inner disk strongly suggest that the starburst in 
NGC 4945 is ongoing and that a post-starburst stage 
of evolution is not yet reached. 
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Fig. 1. (a) 1.4 GHz continuum superimposed on a UK Schmidt Telescope SRC (Science Research Council; now SERC 
= Science and Engineering Research Council) J (yellow sensitive emulsion/filter combination) image. The contour 
levels have flux densities of 5, 7.5, 10, 15, 20, 25, 30, 40, 60, 120, 250, 500, 1000, 2000, 3000, 4200 mJy beam" 1 . The 
resolution is 19" (R.A.) x 25" (Dec.) (i.e. an average of ^23") along major and minor axes, (b) Integrated Hi emission 
superimposed on an SRC-J image. The contours denote —5000, —500, 3, 7.5, 15, 30, 50, 70, 90, 99% of the integrated 
emission peak of 3.8 Jykms -1 beam -1 . The resolution is as in (a), (c) Distribution of CO(2— 1) emission integrated 
over velocity superimposed on a SRC-J image. The contours are 10, 20, 40, 80, 160, 320, 640, 1280 K km s -1 beam -1 
on a T m b scale. The resolution is ~23", and the rms noise is ~ 3Kkms -1 beam -1 . The four white crosses served as 
reference positions for the overlay. 
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Fig. 7. Position- velocity diagrams of Hi with a resolution of ~23" x 7.25 km s -1 , obtained along the major axis at 
(a) PA — 49° and (b) PA = 43°. The velocity axis is labeled with respect to 560 km s" 1 . The contour levels have 
flux densities of -500, -100, -25, 4, 7, 10, 15, 20, plus (a) 30, 40, 50, 55, and (b) 30, 40, 45 Jy beam- 1 . 
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Fig. 8. (a) H I velocity field obtained from gaussian fitting to observed spectra. The iso- velocity contours are at 
385, 410, 435, ... 710, and 725 kms -1 , (b) Brandt model velocity field for Hi. The contours are spaced at 25kms~ 
intervals, (c) Residual velocities, after subtraction of the model field from the observed field. The contours are spaced 
at 20kms _1 intervals. 



